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ABSTRACT 


Name of student: Himanshu Chaudhary Roll No; 9810512 
Degree for which submitted; M.Tech. Dept. ;Mechanical Engg. 

Name of thesis supervisor: Prof Prashant Kumar 

Month and year of thesis submission: January, 2000 

A combined experimental and finite element method is exploited to determine 
d\ namic interlaminar fracture toughness of glass fabric/epoxy laminates. A 7-8 mm thick 
specimen made of 64 laminae and having a precrack of length 40 to 50 mm at the midplane of 
one of its edges is used in the present study. Tire front face of the specimen is bonded to a 
rigid block and the rear cantilever is screwed to a load bar. In the load bar, a compressive 
clastic stress pulse is generated by impacting it with a striker bar. When the incident 
compressive pulse reaches the specimen, a part of energy is reflected into the load bar and the 
rest passes to the specimen. The incident and reflected stress pulses are monitored by two 
strain gauges on the bar, which determine the deflection history of the cantilever end. The 
encrgv' transferred to the specimen initiates growth of the precrack at high speed under the 
impact conditions. The propagation of crack tip is monitored by two strain gauges of width 
0.9 mm and 2 mm length bonded ahead of the crack tip on the side face of the cantilever end. 
The first strain gauge is bonded very close to the crack tip. The crack velocity determines 
initiation time of crack propagation. 

The experimental results, the deflection history' of the cantilever end, the crack 
\ elocit\'. the initiation time of crack propagation, as w'ell as material properties are used as 
input data in a finite element (FE) code to evaluate d>’namic interlaminar fracture toughness 
of the laminates in terms of J-integral( j ). A model of gradual release of nodal forces 
simulates the propagation of the interlaminar crack. The FE code gives the variation of J- 
integral with time. Two sets of specimens of having different crack lengths are investigated to 
stud>’ the effect of crack length. First set of the specimens employs a short crack of 40 mm 
while the other employs a long one of 50 mm. For both sets of experiments the velocit\' is 
kept between 660 and 2200 m/s. The initiation fracture toughness for the short crack length is 
found to be between 55 and 536 J/m^ and that of for the long crack length between 38 and 400 
J/my On the other hand, the values of propagation fracture toughness lies between 6 and 18 
J/m' for short crack length and between 3 and 20 J/m^ for long crack length. The effect of the 
crack length was found insignificant. 
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Chapter 1 


Introduction 


1.1 Introduction 

Glass, carbon, or kevlar-fiber-reinforced polymers (GFRP, CFRP, and KFRP) 
have emerged out as an important class of engineering materials. Since they offer 
outstanding mechanical properties (strength, stiffness and lightness), unique flexibility 
in design capabilities, and ease of fabrication, they are increasingly used as structural 
members. Additional attractive properties, they offer, are corrosion resistance, 
impact resistance, excellent fatigue strength and environmental stability. Therefore, 
though expensive, are finding increasing use in aerospace, transport and sporting 
goods. In addition, there are many opportunities for their wider application in other 
fields like hiking equipment, medical goods and even apparently insignificant things 
like speculator frames. 

Most FRP composites used in structural applications are in laminate form. The 
laminate is commonly made by stacking prepags over each other and curing them at 
high temperature and pressure. High strength man made fibers (carbon & graphite, 
Kevlar, glass, boron) impart stiffness and strength in the plane of the ply but plies are 
bonded to each other with comparatively low strength matrix material; fiber do not 
reinforce the laminate through thickness. However, interlaminar strength is usually 
adequate for an FRP structure loaded under quasistatic conditions. Glass fibers are the 
most common of all reinforcing fibers for polymer matrix materials because of low 
cost and high strength although they have comparatively low modulus and poor 
abrasive resistance, which decreases its strength and adhesion to polymer matrix resin 
in presence of moisture. 

The strength of an FRP laminate comes into question when a foreign body is 
impacted on it at subsonic velocity. An intense damage is observed at the center of the 
impact with fiber breakage, matrix cracking and fiber pull out. Further, the damage 



spreads laterally to the considerably large area mainly through interlaminar separation 
even if the impacting body is of Low Mass (10-30 g) (Kumar and badri, 1993). On the 
other hand, when a foreign body impacts a metallic sheet, material at the impact 
center yields, resulting in formation of a small dent. Therefore, in such cases the FRP 
laminates are less tough than conventional materials. The speed of expanding fronts 
of delamination cracks is as high as 200-500 m/s (Takeda et al .,1982 ). When glass 
fiber reinforced epoxy laminates are impacted it was found that if the total 
delamination area between various plies is multiplied by the quasistatic energy release 
rate, the net energy exceeds the energy of the impacting mass (Kumar and 
Naraganan, 1993). The gives an idea that the interlaminar cracks at high speeds 
propagate at lower toughness resulting in a fairly large damage area. Therefore, the 
characterization of the interlaminar fracture toughness under dynamic creak 
propagation is of prime concern of FRP laminates. 

Experiments to determine dynamic fracture toughness have been mostly done 
on large plates having a through the thickness crack. This work is generally carried 
out through the method of caustics. Although a considerable amount of work has 
been done to study the fracture phenomena through numerical methods but only large 
rectangular double cantilever beam were given attention. Owen and Shantaram 
(1997) introduced the use of finite element method to study dynamic crack growth. 
Nishioka and Atluri (1982a, 1982b) investigated the crack propagation and arrest in a 
high strength steel DCB specimen using moving singular dynamic finite element 
procedure. 

The study of fast moving delamination cracks in a laminate, whose thickness 
is usually less than 5 -8 mm, is complex from an experimental and numerical point of 
view. Modeling of dynamic propagation requires having measurement within a very 
short period during crack propagation. In case of static case the crack propagates 
when the fracture toughness exceeds a certain critical value. Crack propagation 
criterion for dynamic case is not well defined. Further, studies of dynamic crack 
propagation require accurate loading and measuring techniques. These techniques are 
readily not available with commercially available facilities. Keeping these views in 
mind, a combined technique is developed, which makes use of a numerical simulation 
using experimental data from a less sophisticated instrument. 



1.2 Literature Survey 


A large number of investigations have been carried out on dynamic crack 
propagation employing a through-the-thickness crack in a large plate. Stress waves 
are the continuously generated at the crack tip of fast moving crack. As these stress 
waves reach the boundary of the specimen, they are reflected. Some of these reflected 
stresses waves come back to the region infront of crack tip. Experiments with large 
size plate are controlled such that the monitoring of important parameters is complete 
before the reflected waves arrive. In other words, the experiments are controlled to 
isolate the effect of specimen edges. In case of an interlaminar crack in a fiber 
glass/epoxy laminate, which is usually not thicker than 6 to 8 mm, the free surfaces 
are always close to the crack tip and there is no way one can isolate the region of 
interest from the reflected waves. Therefore, a new approach is required to find 
dynamic toughness of interlaminar cracks. 

An important aspect of dynamic crack propagation study is the measurement 
of crack velocity. Takeda et.al, (1982) used high speed photography to measure the 
interlaminar crack propagating speed in composite laminates. They observed that the 
interlaminar crack in composite laminates moves at 200-500 m/s, when impacted by 
foreign bodies. 

Several groups did initial and extensive investigations on crack propagation of 
through the thickness crack in a large plate. Ravi Chandar and Knauss (1982, 1984a, 
1984b, 1984c, 1984d) have done an exhaustive study of dynamic crack propagation 
phenomena using the method of caustics for Homalite-100. This included study of 
crack initiation and arrest, microstructure aspects, crack branching and interaction of 
stress waves with the crack tip. Rosakis, Duffy and Fruend (1984) performed 
dynamic crack propagation experiments on double cantilever beam specimen using 
wedge loading. Zehnder and Rosakis (1990) used optical method of reflected caustics 
combined with high-speed photography to investigate the dynamic fracture initiation 
and propagation in 4340 steel specimen. 



These investigations, based on caustic have been very effectively used to study 
the fracture phenomena in a large plate with through the thickness crack. The method 
is not likely to be effective for the interlaminar crack in a slender sheet because the 
strength of the interlaminar bond is very weak and it will introduce very little strain in 
the bonded sheet making the size of the caustic negligibly small. 

Ravichandran and Clifton (1989) developed a special technique to study the 
initiation and propagation of crack in the steel under dynamic impact loading. They 
presented a plate impact experiment and an associated finite difference model to study 
the fracture process that occurs in sub-micron loading. A disc containing a 
prefatigued edge crack on the midplane upto half way across the diameter is impacted 
by a thin flyer plate of same material. A compressive pulse propagates through the 
specimen and reflects from the rear surface as a tensile pulse of Ips duration. Using 
the laser interferometer system monitors the motion of the rear surface. 

Berger and Dally (1990a) used a series of strain gauges ahead of the crack tip 
at a certain predetermined location to monitor the strain and crack propagation. 
Berger, Dally and Sanford (1990b) also used strain gauges ahead of the crack tip to 
determine dynamic stress intensity factor associated with a propagating crack. 

Nishioka and Atluri (1983) studied the use of path independent J -integral for 
dynamic crack propagation by the finite element method. Other path integrals were 
also investigated along with J -integral. Numerical results showed that combined use 
of J -integral and the finite element method is a useful tool to obtain the fracture 
parameters such as stress intensity factors and energy release rates. 

Kolednik (1991) presented his theoretical study for physical interpretation of 
the J-a curves for elastic-plastic fracture. He derived the difference between the exact 
initiation toughness and crack growth toughness using energy balance under quasi- 
static conditions. He considered three point bend specimens each of that consists of 
two parts glued together along the ligament. The analysis was made on large 
specimens. 



Verma (1995) developed a combined experimental and numerical technique to 
determine interlaminar dynamic fracture toughness under impact loading. An 
interlaminar crack is propagated at very high speed in DCB specimen made of two 
steel plates which are bonded together by epoxy with a precrack. Strain gauges of 
gauge length 0.2 mm are bonded ahead of crack tip for monitoring crack propagation 
speed. Using the experimentally obtained data (deflection of steel plate, initiation 
time and crack propagation history with time) to FE code J-integral is determined. 
Initiation toughness and propagation toughness are evaluated. 

Sun and Grandy (1988) investigated dynamic delamination fracture toughness 
in a [90/0]5 s T900/934 graphite/epoxy laminate using impact loading. Delamination 
cracks of different sizes were embedded at the midplane of the composite specimen. 
The threshold impact velocity that causes propagation of delamination crack was used 
in the dynamic analysis with the finite element method. From the finite element 
solution, the time-history of the strain energy release rate was calculated. The critical 
strain energy release rate was taken equal to that of maximum value of the response 
history. 


Raman P. Singh et al. (1996) described various experimental observations for 
dynamic intersonic decohesion of bimaterial interfaces. Two separate but 
complementary optical methods are used in conjunction with high-speed photography 
to explore the nature of the large-scale contact and shock wave formation at the 
vicinity of running cracks in two different bimaterial systems. 

R.W. Truss et al. (1997) had attempted to find interlaminar and intralaminar 
fracture toughness of uniaxial continuous and discontinuous carbon fiber/epoxy 
composites, using compact tension and double cantilever beam test geometry. The 
discontinuous carbon fiber/epoxy composites have been found to have a slightly 
misalignment of the fibers from the average fiber direction and this misalignment was 
found to increase both initiation fracture toughness and to greater extent the 
propagation fracture toughness. The increase in fracture toughness in discontinuous 
carbon fibers/epoxy samples was due to fibers bridging the crack and this has been 
modeled as if the the bridging fibers provide an increase in compressive stress across 



the crack analogues to a craze at the crack tip. 


John Lambros and Rosakis A.J (1997a) investigated dynamic delamination of 
thick fiber reinforced polymeric matrix composite laminates using optical techniques 
and high-speed photography. They had used 65% fiber volume fraction of 
Graphite/epoxy laminates consisting of 48 plies. Square plates ofl52 mm X 152 mm 
dimensions were impacted in an out-of-plane configuration using high speed (Im/s - 
3 m/s) gas gun. Real time imaging of the laminate out-of-plane displacement was 
performed using lateral shearing of interferometer of Coherent Gradient Sensing in 
conjunction with high speed photography. Delamination speeds up to 1800 m/s were 
observed. 

John Lambros and Rosakis A.J (1997b) investigated dynamic crack initiation 
and grov\4h in unidirectional graphite/epoxy plates with high-speed photography. 
Edge notch plates are impact loaded in one point bend configuration using drop 
weight tower. Initiation fracture toughness data were reported and significant 
dynamic effects are observed through emission of stress waves from the propagating 
crack tip. 

Ramakrishna (1996), Babu (1998) and Mallikharjuna (1998) attempted to 
investigate initiation and propagation toughness of unidirectional GFRP laminates 
under impact loading. In these studies, Verma’s(1995) experimental technique was 
exploited with some changes. Ramakrishna (1996) found that the initiation toughness 
lies between 50 and 250 J/m^ and propagation toughness between 20 and 50 J/m^ 
Babu and Mallikhajuna (1998) used superior specimen and conducted experiments at 
crack velocity varying in the range of 517 to 1000 m/s. They used velocity strain 
gauges of smaller length (0.2 mm). They reported J-initiation between 245 and 620 
J/m^ and J-propagation between 45 and 83 J/m^ Ramakrishna, Babu and 
Mallikharjuna used GFRP laminates to investigate interlaminar dynamic fracture 
toughness. 



1.3 Present Work 


To determine interlaminar dynamic toughness of glass fabric/epoxy laminates 
under impact loading, a technique developed by Verma (1995) are exploited. The 
technique is initially developed on a double cantilever beam specimen made from 
hardened steel strips bonded together with epoxy. One of the cantilevers is impacted. 
A FE code is developed to simulate the wave propagation in the DCB specimen and 
to determine the initiation and propagation toughness of the interlaminar crack. This 
hybrid technique is extended in the present work to determine initiation and 
propagation toughness of specimen made of glass fabric /epoxy laminate. The 
technique employs a combined experimental and numerical analysis. 

The technique used for glass fabric/epoxy laminates consists of applying an 
impact load on a specimen with a precrack (Fig. 1.1). The portion that is on front side 
of the precrack remains straight during experimentation because it is bonded to a solid 
block. The rear portion of the specimen acts as a cantilever plate. Strain gauges are 
mounted ahead of the crack tip to monitor the crack velocity. The experimental work 
determines the end deflection of the cantilever, crack velocity and crack initiation 
time. These are used as inputs to the FE code to simulate the crack propagation. 

Chapter 2 discusses the specimen preparation and experimental setup. The 
brief details of finite element code are covered in Chapter 3. Results of the present 
work and comparison of the same with that of other researchers are presented in 
Chapter 4. Chapter 5 concludes the present work and scope for further work. 




Fig. 1 . 1 Specimen under impact loading 


Chapter 2 


Specimen and Experimental Set up 


In the present study, a combination of experimental and a finite element 
package is employed to find interlaminar dynamic fracture toughness of polymer 
composites. Schematic diagram of the experimental setup is shown in Fig. 2. 1 and 
the photograph of the setup in Fig. 2.2. Front face of the specimen is bonded to a 
rigid support and the rear cantilever is screwed to the load bar. The load bar is 
impacted with a striker bar. An air gun using compressed nitrogen gas is used for 
acceleration of the striker. Strain gauges mounted on load bar at an appropriate 

location measure the stress pulses propagating in the load bar. To monitor crack 
propagation two velocity strain gauges are bonded on side fkce of the DCB specimen. 

In this section specimen preparation, experimental setup, experimental 
measurements and data processing are discussed. 

2.1 Specimen 

The specimen used in the present work is glassfibre fabric reinforced epoxy 
composite. This section describes the specimen geometry, raw materials, and 

preparation of laminates and preparation of specimen from laminate. 

2.1.1 Specimen Geometry 

The geometry of the specimen is shown in Fig. 2.3. Figure 2.4 shows the 
orthographic views of the same. It is 25 mm wide and 180 - 200 mm long with an 
edge precrack on mid plane of the specimen. The precrack is introduced during the 
fabrication of the specimen. The portion of the specimen, which is on front side of 
the precrack is cut and removed so that the load bar can impact the cantilever . A 
small hole of diameter 6.2 mm is drilled on the cantilever so that it can be screwed to 
the load bar with an M 6 cap screw and a specially designed washer. 


2.1.2 Raw Materials 


Fiberglass fabric and epoxy resin mixtures are used as raw materials of the 
specimen. Composition of epoxy resin mixture is as follows : 

Araldite LY556 1 00 parts by weight 

Hardener HT976 35 parts by weight 

Accelerator XY73 ] part by weight 

Coupling Agent 0.5 part by weight 

(y- Amino prophyl triethoxy silane) 

Fiberglass fabric, purchased from M/s HARSH-DEEP INDUSTRIES, 
Ahemdabad, is used in the present work. The properties of fiberglass fabric are 
dependent on fabric construction, i.e., the number of yam per cm in each direction, 
weave pattern, and yam type. Specifications of the fiberglass fabric used as 
reinforcement material in the present work are as follows: 

Number of yam per 25 mm in warp = 40 
Number of yam per 25 mm in weft =36 
Area density =146 gm/m^. 

2.1.3 Laminate Preparation 

Prepregs of fabric and epoxy resin is made through a 150 mm wide prepreg- 
making machine available in the laboratory. It is prepared at 95°C and stored at - 
18°C. For making a laminate, this prepregs is cut for required approximate length of 
240 mm. Each cut lamina is stacked together to get required thickness; 60-64 lamina 
are generally required for a specimen of 7 to 8 mm thickness. Precrack length, 35 to 
65 mm, is introduced by inserting a thin teflon film of 20 pm thickness at the 
midplane of the laminate at one of its edge. The stacked laminate is cured in 
hydraulic press at 120° C temperature. Initially laminate is pressed between two 
plates at gradually increasing pressure, 23.7 kPa per minute for half an hour. Then 
pressure of 710 kPa and temperature 120° C are maintained for 1 hour. Finally, 
maintaining the same pressure and increasing temperature to 150oC for two hours 
does post curing. The pressed laminate is then allowed to cool to room temperature. 
A diamond cutter cuts this cured laminate in required width of 25 mm. 



2.1.4 Specimen from Laminate 


Now the specimen is marked for drilling a hole to have the desired crack 
length. After drilling the hole of 6.2 mm diameter, the front portion of the specimen 
near the hole is removed through a saw. Furthermore, the front face of the specimen 
is serrated with a hacksaw for better bonding to the solid block. Final shape of the 
specimen is shown in Figs.2.3 and 2.4. 

In the present study the specimen used is 7-8 mm thick, 25 mm wide and 220 
mm long whereas crack length for one set of experiments is closed to 40 mm and for 
other set of experiments approximately close to 50 mm. The portion on the front side 
of the precrack remains straight during experimentation because it is bonded to the 
solid block. The rear portion of the specimen works as a cantilever. 

2.2 Experimental Setup 

Main integral components of experimental setup include air gun, striker, load 
bar, bridge circuits, digital storage oscilloscope and solid block (Figs. 2. land 2.2). 
The slotted face of the specimen is bonded with araldite to the rigid block with the 
help of C clamps and the block is fastened to base of set up with two Ml 2 bolts and 
washers. The load bar is screwed to the rear end of the specimen with a M6 screw 
lUa washer. V'ihen the load bar is impacted by a striker bar, from air gun, elastic 
compressive and reflected tensile stress pulses are generated in it. These elastic 
stress pulses are monitored by two strain gauges bonded on the bar in the longitudinal 
direction. They are connected to an oscilloscope through a bridge circuit. The two 
velocity strain gauges bonded ahead of crack tip are also connected to the 
oscilloscope each through another bridge circuit. The complete bridge circuit is 
discussed subsequently. 

2.2.1 Striker and load bar 

A compressive stress pulse is applied to the specimen through load bar, striker 
and air gun of nitrogen. The load bar is properly aligned with the centerline of barrel 
of the air gun. The striker bar is accelerated in the air gun and its impacting face is 



spherical so that impact would be at center of end face of the load bar. Striker and 
load bar are made of cold rolled mild steel of 19 mm diameter. V.^en the striker, 
from air gun, hits the load bar a longitudinal compressive stress pulse is generated in 
the load bar. As this compressive stress pulses reaches the specimen end of the load 
bar, part of it transmitted to the specimen and remaining is reflected back as tensile 
pulse. 


Two strain gauges bonded at diametrically opposite locations on the load bar 
pick up these stress pulses. The reponses of the gauges are recorded in a storage 
oscilloscope through a bridge circuit. 

2.2.2 Stress Pulses 

V.'ihen the striker bar impacts the load bar, the compressive incident pulse (cti) 
propagates towards the specimen and is recorded at location 1 as shown in Fig.2.1. 
The reflected tensile pulse (as) is recorded at the same location but at different time 
and is denoted by point 3 in the t - x diagram (Fig. 2.5). A typical record of the 
incident and reflected pulses in load bar is shown in Fig.2.6. 

To determine displacement of cantilever end one-dimensional wave 
propagation theory for elastic wave propagation in the load bar and the striker bar is 
used. /V.^jave propagation in load bar is analyzed through relation along characteristic. 
Line 1 - 2 is the characteristic along positive direction and relation between the stress 
and particle velocity is given by 

da - pcdv = 0 (along the positive characteristic, 1-2) 

dp + pcdv = 0 (along the negative characteristic, 2-3) 
where a is stress, v is particle velocity and p is density of the load bar material, c is 
the longitudinal wave velocity in material of the load bar. By using above equations 
the relation along characteristic 1-2 simplifies to 


02 - PCV 2 = oi - pcvi. 


( 2 . 1 ) 


Along the characteristic 1-5 


Oi + pcvi = 05 - pcv5. (2.2) 

However, oj = 0 and V5 = 0 because the load bar is initially at rest and stress wave 
never reaches at point 5. The above equation gives 


a] = - pcv]. 


(2.3) 


Substituting Eq. 2.3 in Eq. 2 . 1 , one obtains 

02 - pcv2 = 2ai. (2.4) 

Relation along characteristic 2-3 and 3 - 4 are 

03 + pcV3 =02 + pcv2 (2.5) 

03 - PCV3 = 04 - PCV 4 . (2.6) 

However, 04 = 0 and V4 = 0 at point 4 because the striker and load bar are of same 
material and diameter and the striker comes at rest. The above two Eqs. Yield 


O2 + PCV2 = 203. 


(2.7) 


O, -a, 


pc 


(2.8 


It is to be noted that nature of Oj is compressive and 03 is tensile and therefore 
the particle velocity of the cantilever end can be found by taking sum of the absolute 
value of incident and reflected pulse. Displacement U2 at point 2 with time t is found 
by integrating Eq. 2.8 


pc 


(2.9) 



The time, t, at cantilever end is measured only after the head of the incident pulse 
reaches the strain gauge of the load bar. Similarly time t, is measured for reflected 
pulse when head of the reflected tensile pulse reaches at the strain gauge of the load 
bar. 

2.2.3 Details of Bridge Circuit 

The bridge circuit converts the change in resistance of the strain gauges 
bonded to the load bar into potential difference. This potential difference is recorded 
by oscilloscope. Figure 2.7 shows the configuration of the bridge circuit used in the 
present study. Rj, R2, R3 and R4 are four strain gauges each having a resistance of 
120 ± 0.3 £1. Riand R 2 are mounted on the load bar whereas R 2 and R 4 represent 
dummy gauges. The dummy gauges are also mounted on a 19 mm diameter steel bar. 
The circuit is balance by connecting a 1.0 fl resistance in series and two variable 
resistors 100 kD and 1.0 KQ parallel to one of the dummy gauges R2 and R4. 
Balancing of bridge is necessitated to give zero voltage output for no change in 
resistance of the active strain gauges. A calibration resistance Rc is connected parallel 
to one of active gauge (Ri in diagram) through a switch K. The value of calibration 
resistance is 47 kfl. For a balance bridge circuit the output voltage (Ae) is given by 


R,R,E 

(R,+R,)= 


AR^ 

R, 


AR; AR3 AR, 

R3 Rj R4 


where E is the voltage applied (9 V in this study) to the bridge circuit. 


(2.10) 


This equation shows that similar (both positive or both negative) changes in 
resistance of opposite arms of the bridge circuit are added up and dissimilar (one 
negative and other positive) changes are cancelled out. Thus, active strain gauges at 
opposite arms only records compressive pulse and bending effect is neglected. 


The relation between the strain in strain gauge and corresponding change in its 
resistance is governing by the following equation 



8 



AR/R 

(GF) 


( 2 . 11 ) 


where GF is the gauge factor and s is the strain recorded by the strain gauge. 
1 he change in resistance of arm CD after connecting R<; is given by 


AR3 


R, 




R- +R 


c J 


( 2 . 12 ) 


which gives 


AR, 


R, 


(2.13) 

R3 R,+R, 

Corresponding to this change in resistance a voltage difference (Vc) will occur 
betw'een terminals A and C. Similar resistance change also occurs in R 3 when the 
load bar experiences the load pulses. Thus the simulated strain (se) due to shunting 
the calibration resistance is given by 


e = 


R, 


2(R3+RJGF 


(2.14) 


The strain corresponding to calibration voltage (Vc). It is worth noting that a factor V2 
in Eq.2.14 is applied to average out the strain recorded by the two strain gauges Ri 
and R?, which are at opposite arms of the bridge circuit. The strain in the load bar 
corresponding to a voltage V recorded on oscilloscope is given as 

E = -^V (2.15) 

Vc 

2.2.4 Digital Storage Oscilloscope 

A 4-channel digital storage oscilloscope of DS01624, Gould Inc., UK is used 
to monitor strain gauges response. It stores data in 50-memory block and can be 
recalled. It has 12 bit resolution is 4096 point on the screen in horizontal and 8 bit 
resolution i.e. 256 points in vertical direction. Therefore per channel storage points 



available in horizontal direction are 1024 and that in vertical direction 64. The 
maximum sensitivity in the vertical direction is 0.025 V and that in horizontal 
direction is 0.25 |is. For 4 channel horizontal sensitivity become 1 [is. The stored 
data are transferred to a personal computer by the GPIB software. 

2.2.5 Rigid Block 

A rigid block of mild steel of dimensions 75mm X 75 mm X 150 mm is used 
to bond serrated face of the specimen. The Serrated face of block is kept towards the 
serrated face of specimen. Two holes in the block are oversize (18 mm diameter) so 
as to adjust the specimen normal to the load bar. The rigid block is bolted to the rigid 
base plate after aligning the specimen to the load bar. Care is to be taken to ensure 
that the crack plane is normal to the load bar; otherwise specimen will be preloaded 
and the numerical simulation would be inaccurate. To use same block repeatedly the 
block is cleaned by burning araldite and then cleaning with acetone. 

2.2.6 Velocity Strain Gauges 

The two type of gauges used are supplied by Tokyo Sokki Kenkyijo Co., Ltd. 
Japan of following specification : 


(I) BFLA-2-8 Resistance : 

120 ± 0.3 n 

Gauge factor ; 

2.11 ±1 

Length: 

2 mm 

Width 

1 mm 

(II) BFLA-5-8 Resistance : 

120 ±0.3 0 

Gauge factor ; 

2.10+ 1% 

Length; 

5 mm 

Width 

1 .3 mm 


Verma et.al (1995), Babu and Mallikharjuna (1998) used the strain gauge of 
type FLG-02-1 1, gauge factor 2.05 and gauge length 0.2 mm. Since Verma (1995) 
used specimen of isotropic steel plate as bonding face, they worked properly. In the 



present study, the material used is heterogeneous therefore 0.2 mm length is unable to 
sense the singular strain peak properly. Furthermore, BELA strain gauges are more 
suitable for this study because they are specially designed for composites. 

The backing film of the gauges is large in area and therefore a portion of the 
gauge is chopped (Figs. 2.8 and 2.9). The gauge is bonded on the side face of the 
specimen such that the top portion of the active portion is 0.5 to 0.8 mm away from 
the crack plane. 

Before bonding the strain gauges, the faces of the specimen to be bonded are 
polished with the help of emery paper. Then it is cleaned by acetone. The first strain 
gauge is placed very close to the crack tip approximately 1-2 mm and the second 
strain gauge 3-5 mm from first the gauge. These strain gauges are perpendicular to 
crack plane. Fewi-quick supplied by Pedilite Industries is used to bond the gauges. 
Position of the gauges with crack tip and plane are shown in Fig. 2. 10. 


2.2,7 Crack Velocity 

Crack velocity is obtained by monitoring the responses of the strain gauges 
bonded ahead of the crack tip (Verma 1995). Verma used to bond three strain gauges 
at a certain angle. In present work, a modification of the above technique is used. 
Only two strain gauges are used and both the strain gauges are kept at an angle 90° to 
the crack plane. It is noteworthy that strain gauges are kept at a short distance (0.5 to 
0.8 mm) from the midplane or crack plane of the specimen. Since the strain gauge are 
very close the crack plane, the singular strain field of the crack tip gives a strain peak 
when the crack tip passes close to the strain gauge. There are two such velocity strain 
gauges; the distance between the two strain gauges is measured accurately through 
travelling microscope; the velocity can be found using formula (Fig.2. 11) 


Velocity = 



t-, -t, 



2.2.8 Initiation Time 


The initiation time could be found out by extrapolating the data procured by 
the crack velocity strain gauges as shown in Fig. 2.11. Extrapolation procedure to 
find initiation time is justified on the following observation. 

• The first velocity strain gauge is very close to the tip of the precrack (in 
the range of 1 to 2 mm). The initiation time of propagation of these kinds 
of experiment is very small. This can be justified from the experimental 
work of John Lambros and Ares J. Rosakis (1997a) who conducted on 
impact loading test on a one-point bend edge-notched unidirectional 
graphite/epoxy composite plates. For obtaining crack initiation and 
monitoring crack growth, on experimental technique with high speed 
photography was used. The measured crack velocity Fig. 2. 12 is in 
microsecond domain; the crack, after initiation, acquires high velocity 
immediately with rise time much less that one ps. The experimental 
results show that the extrapolation invoked in this study in finding 
initiation time work reasonably good. It is observed clear from Fig.2.12 
that crack velocity does not change drastically after initiation time. 

• The crack velocity of interlaminar crack is very high (600 to 2000m/s) and 
it takes a short time for the crack to move to the first velocity strain gauge. 

2.3 Experiment Procedure 

The specimen bonded with velocity strain gauges is screwed to the load bar 
through M6 cap screw and the front plate is bonded to rigid block by araldite. To cure 
bonding materials, block and specimen together is kept for 12 hours under bulb 
heating. Leads of velocity strain gauge are soldered to bridge circuit arm as discussed 
in section of bridge circuit; bridge circuits are then balanced. Experiments are 

carried out as follows: 

Approximately 3.5 bar pressurized gas is filled in air gun from cylinder. 
Before releasing pressure it is made sure that oscilloscope is checked for normal and 
storage mode. The pressurized air is released to accelerate the striker, which impact 
the load bar. A compressive stress induced in the load bar propagates through strain 



gauges and reaches the cantilever of the specimen. A part of the compressive pulse is 
transmitted to specimen that propagates the crack. Remaining pulse is reflected back 
as tensile pulse that again is recorded through the gauges of the load bar. The 
oscilloscope records strain pulses in the load bar, which are analyzed using one 
dimensional elastic wave analysis to obtain the displacement history of the cantilever 
end. Responses of the strain gauges on cantilever face, recorded by the same 
oscilloscope in other channels, yield the initiation time of the crack growth and crack 
velocity. 

This whole set up was initially developed by Verma at.al (1995) for 
investigating fracture toughness of bonded steel plates. It has now been modified to 
determine interlaminar toughness of fabric glass fibers reinforced epoxy laminate 
under impact loading conditions. 

2.4 Closure 

The experimental set up measures the parameters: 

1 . Displacement vs. time relation of the cantilever end 

2. Interlaminar crack velocity and 

3. Initiation time of crack growth. 

The experimental data are supplied to the finite element simulation to 
determine toughness required to initiate the crack (Jini) and toughness of propagating 
crack (Jprop).The details of FEM simulation are discussed in Chapters. 



SPECIMEN 



Fig. 2.1 Schematic diagram of the Impact loading setup. 










Fig.2.4 Orthographic views of the specimen 









Fig. 2. 7 Configuration of bridge circuit 




Fig. 2, 10 Locations of velocity strain gauges with reference to the crack plane 





Chapter 3 

Numerical Analysis 


The experimental measurements gives the profile of the displacement applied 
at the cantilever end, the time initiation of crack propagation and crack propagation 
history. These data are used in a finite element code to simulate the dynamic fracture 
behavior of the DCB specimen. Dynamic fracture phenomena has several important 
features. As the crack extends and the boundaries of the body change with time, 
mathematical models are more complex than those of the static case. From 
experimental point of view, parameters are transient and they have to be measured at 
high speed accurately. 

Keeping all these points in mind, Verma at.al (1995) developed a FE code for 
finding J-integral variation with time for interlaminar growth crack between two 
bonded thin plates made of isotropic steel. This FE code is modified by Babu (1998) 
for evaluating initiation and propagation toughness of unidirectional FRP laminates. 
Further the program is modified, in the present study, for fabric reinforced polymer 
composite material. 

This chapter briefly describes the finite element formulation, path independent 
integral, the crack opening scheme and the flow chart of the code. 

3.1 Finite Element Formulation 

3.1.1 Formulation 

The governing system of finite element equations for the linear dynamic 
response of an elastic body is given in the matrix form as 


[M]{U} + [K](U}= {R} 


(3.1) 



where [M] and [K] are the mass and stiffness matrices respectively and {R} is the 
external load vector, and {U}, {U} are the displacement and acceleration vectors of 
the finite element assemblage respectively. Damping forces have been neglected in 
the formulation. In the above equation global matrices may be obtained as 

[M] = Z[M]" 

[K] = Z[K]" 

[R] = S[R]^ 


where element mass, stiffness and traction matrices are given by 

[M]" = Iv p [Mf [N] dV 
[K]^ = Iv [B]^ [D] [B] dV 
{Rr = L [N][T]dS. 

In these equations [N] is shape function, [B] derivative of shape functions, [D] elastic 
constitutive relation matrices and [T] traction vector. 

3.1.2 Mode Superposition Method 

Equation 3.1 can be solved either by time integration or by mode 
superposition; the former is preferred for simulating wave propagation problem in the 
specimen. In this integration scheme there are many different methods, which can be 
classified as “explicit” or “implicit”. Relative advantages and disadvantages of each 
method are given in Bathe (1990). Verma at al. (1997) developed a FE code for 2-D 
finite element discretization in space analysis using the Newmark integration method 
for the time variable (Bathe, 1990). Since width of the specimen is large in 
comparison to the thickness of the specimen, this is a plane strain problem. 

Specimen front face that is bonded to a rigid block (Fig. 1.1) remains straight. 
It is therefore expected that effect of this portion on J-integral is negligible in fact. 



Verma (1995) simulated FE analysis in the entire specimen and showed that the effect 
was insignificant. Later, Mallikharjuna (1998) also verified this. Consequently, only 
cantilever plate of the specimen is analysed to avoid unnecessary computations. 

3.1.3 J-Integral 

Nishioka and Atluri (1983) studied the use of path independent J-integral 
using finite element method for dynamic crack propagation. In this present work, ji 
given by Kishimoto, Aoki and Sakata (1980) is used which is defined as (Fig. 3.1) 


+ PU,u,,dV. 


where W is the strain energy density, Ui is the displacement, ni is the direction cosine 
of the unit outward normal, T; is the traction. Pc', F, Fc"^ aire defined paths and V.Vg 
are volumes. 


J-integral is independent of the choice of the path F (Fig. 3.1) under steady 
state crack growth conditions. In many cases, though these conditions are obviously 
not met fully, the path independence of J -integral can still be established within 
certain allowable errors. To minimize errors, the integral is evaluated over a path far 
away from the crack tip. Then the analysis is relatively insensitive to the finer details 
of the crack tip region. The path is held stationary as crack tip extends in a self- 
similar manner. 

3.1.4 Constitutive Relations for Composites 

From the mechanics point of view fibre composites are among the class of 
material called orthotropic materials. Generalized Hook’s Law for two-dimensional 
orthotropic laminate is given by (Agarwal, B.D., 1990) 

rtj = Qij 8j 

where CTi are the stress components, Qjj is the stiffness matrix, Sj are the engineering 

strain components and i.j = 1,2,.... ,6. 
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Constitutive relations for polymer composites are expressed for three-dimensional 
cases with nine independent elastic constants. However, only four elastic constants 
are needed for two-dimensional orthotropy (Agarwal, B.D., 1990). 


Consider an orthotropic lamina (Fig. 3.2) with the reference axes coinciding 
with axes of symmetry and designed as the longitudinal direction, L and the 
transverse direction, T. Engineering constants for this lamina are elastic module in 
longitudinal and transverse directions, El and Et respectively, the shearing modulus 
associated with the axes of symmetry, Glt, the major Poisson ration vlt, giving 
transverse strain caused by a longitudinal stress, and the minor Poisson vtl, giving the 
longitudinal strain resulting from a transverse stress. 


Q- = 




Q22 




^LT^T 

1- Vlt Vtl 


Relation between the five engineering constants and four independent elastic 
modulus are expressed as follows (Agarwal, B.D., 1990). 


Qn = 


El 


1 


Qi6 - Glt 


In these equations, a relationship exist between four of the five constants 



Vtl 




In the present study fibre glass fabric is used as reinforcift^%^efiaP‘(tlg.*,S.3').-‘ /Since . e 
number of yarns per cm in each principal direction is a^pr 6 '?^ihiStel 5 ^'Jsa^^ev cne;c.an^^^^^ 
assume 

El = Et 
Vlt =Vtl 

Modulus of elasticity and Poisson ratio are directly obtained’ ifeiW tensflec tests s 
conducted in Experimental Stress Analysis Laboratory and' AfcMS 6'# OT Kanji^irr 
The Glt is obtained from expression ( Agarwal, B.D., 1996). 
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As elastic modulus in longitudinal and transfers (El & Et) are erjiial/ffiepefore Glt^s 
given by 

1 1 f 2(l-v^r) ; 1 ] 

E* 4| E, G„j 


where E 45 is elastic modulus for a 45° ply laminate. 

Experimental values for the fibre glass fabric composite obtained from 

experiments are as follows (Appendix A): 

El = Et = 26.00 GPa 
Vlt “ vtl“0.21 
Glt “ 3.5 Gpa . 

These values are used in all experiments of the present study. 


3.2 The FE Code 

3.2.1 Flowchart 

The flow chart of the finite element code developed by Verma et.4 (1995) is 
shown in Fig. 3.3. This programs evaluates J-integral with time. for specified input of 



material properties (engineering constants), crack propagation history, cantilever end 
deflection, mesh size, time step and the path of J-integral. 


3.2.2 Mesh Size 

Seron et.al (1990) showed that higher order elements are not required in the 
analysis of determining J-integral of dynamic problems. Four noded isoparametric 
elements are used in the present study. To take into the consideration of bending as 
well, fine mesh is used. To avoid nonuniformity of mass distribution, which leads to 
reflection of waves, uniform mesh is used through out the specimen. In the present 
work 0. IhxO. 1 125h mm mesh is used where h is thickness of cantilever (Fig.3.4). 

3.2.3 Crack Opening Scheme 

In order to simulate the crack propagation, a gradual nodal release method was 
used. The actual crack (Fig. 3.5) is located at C between the nodes B, D and the 
holding force F at node B is gradually reduced as the crack is advanced to point D. In 
most of the studies reported in literature the force F reduces to zero when the crack 
reaches the next node D. However, it was found that at very high crack speeds, close 
to the Rayleigh wave speeds, a better model is required. For obtaining a gradual and 
smooth opening of crack, Verma et.al (1995) employed an alternative method. The 
holding back force F at the node B is linearly decreased to zero when the crack 
reaches the end of the next element, point E, in Fig. 3. 5. When the crack tip is in 
between nodes B and D, the holding force Fb at node B given by 



where Fbc is the force at node B when the crack tip was at point B, bi is the crack 
extension beyond node B and d is the element length as shown in Fig. 3.5. When 
crack propagates beyond the node D to point Di, the holding force at node B is given 
by 



3.2.4 Time Step 


_ I d +b; 
^Bc _ 2d 


Time step influences the analysis if not chosen properly. The choice of time 
step (At) for time integration is important for an accurate solution. An optimum 
choice of time step is At = d/c, where d is smallest mesh size and c is the fastest wave 
velocity (Verma et.al., 1995) 

It is observed that upto initiation time, J-integral does not vary much on time 
step. Fig. 3.7 shows the variation of J-integral for three different time steps, 0.2, 0.05, 
0.02 ps. It is observed that upto initiation time, J-integral does not vary much on 
time step. Fig.3.7 shows that amplitude of oscillations in propagating phase varies 
with time step. The amplitude is more for higher value of time step. Therefore a 
coarse time step of 0.2 ps is selected prior to the growth of the crack to reduce number 
of iteration and fine time of step of 0.02 ps is used after the crack growth is initiated. 
By choosing two time steps, overall iterations number used for running program are 
reduced reasonably 

3.2.5 Initiation and Propagation Toughness 

Impact energy, imparted at the end of the rear cantilever, propagates within the 
specimen and reaches the crack tip. J-integral starts building up to initiation time and 
the crack does not advance (stationary crack). Fig.3.8 shows the variation of J-integral 
with time for stationary crack. J-integral builds up sharply, attain a peak and then 
decrease to a low value. However, a stage reaches when the crack tip cannot 
withstand the loading at the crack tip and the crack starts propagating. The value of J- 
integral at which the crack starts moving is known as initiation toughness (Jini). As 
soon as the crack starts propagating J is reduced drastically to a much lower value and 
after some time, it stabilizes (Fig. 3.9). This stabilized value of the J-integral is 
treated as the propagation toughness (Jprop)- 



3.3 Closure 


The aim of this Chapter was to describe the finite element formulation, the 
path independent integral, the crack opening scheme, the initiation and propagation 
toughness and the summary of the program. 

Initiation and propagation toughness are determined through the technique that 
combines experimental results with numerical simulation, discussed and tabulated at 
length in chapter 4, 
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Fig. 3.3 Glass fabric 



Fig. 3.4 Flow chart of the FE code 
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Fig. 3. 8 Atypical variation of J-integral for stationary crack 
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Fig. 3.9 A Typical variation of J-integral for propagating crack 




Chapter 4 


Results and Discussions 

4.1 Introduction 

Experiments are conducted to determine the deflection of the cantilever end, 
crack velocity and initiation time. These are the inputs along with the material 
properties to the FE code as discussed in Chapter 3. The numerical simulation 
determines variation of J-integral with time for stationary and propagating crack. The 
chapter briefly describes the results of initiation toughness, propagation toughness and 
their comparison with the dynamic toughness of interlaminar crack investigated by 
others. 

4.2 Details Procedure of Analysis 

In this section, a detail analysis of Experiment S-1 is discussed. When the 
specimen is impacted by the load bar, a part of the energy is transmitted to the 
specimen. This energy deflects the cantilever. The remaining energy of the load bar 
returns back as a compressive pulse. Strain gauges bonded on the load bar monitor 
the stress pulses. The response of the load bar strain gauges and the velocity strain 
gauges bonded on the side face of cantilever as recorded in the oscilloscope are shown 
in Fig. 4.1(a). A zoomed view in horizontal direction of the oscilloscope records is 
shown in Fig 4.1(b). 

The velocity of the load bar end V 2 is evaluated by one dimensional wave 
propagation theory as discussed in Sec.2.2.2 ^ven by the formula. 


where 01 = compressive incident pulse 
03 = tensile reflected pulse 
p = density of material 
c = longitudinal wave velocity 



The variation of velocity with time is shown in Fig. 4 . 2 . The end displacement U 2 of 
the load bar is obtained by integrating the velocity history (Sec.2.2.2) using 



^ 3-^1 

pc 



Displacement of the load-bar-end is equal to the deflection of the cantilever 
end because load bar is screwed to the specimen. The variation of deflection of the 
cantilever end with time is shown in Fig. 4.3. The loading of cantilever starts at the 
time corresponding to midpoint of the head of the incident pulse and head of the 
reflected pulse. This is denoted as REFERENCE in the Fig. 4. 1(b). 

When the crack tip passes the velocity strain gauges, the singular strain field 
gives a strain peak. The magnified view of the response of the velocity strain gauges 
is shown in Fig. 4.4. The time difference between the peaks of strains and the distance 
between location of the gauges give average velocity between strain gauges. The 
crack length as a function of time is plotted in Fig.4.5. The initiation time is obtained 
by extrapolating the crack velocity upto precrack length. 

The experimental data (the deflection of cantilever end, crack velocity and 
initiation time) are used as input to the FE code as discussed in Chapter 3. The 
variation of J-integral in stationary and propagating phases are shown in Figs.4.6 and 
4.7. The rise of the interlaminar toughness until the initiation is also obtained 
experimentally by Lambros John & Rosakis A.J. (1997). They investigated the 
variation of dynamic fracture toughness upto the initiation time on unidirectional 
Graphite/epoxy composite plates with edge notched under the impact loading (Fig. 
4.8). 

Once the FE code is executed upto the initiation time, the crack propagation 
module is called (Verma, 1995). The location of the crack tip in each iteration is 
known through the experiment; an appropriate factor, as discussed in Sec.3.2.3, is 
applied to the holding back force for modeling crack propagation. As discussed in 
Sec. 3.2.4 a time step of 0.02 ps is chosen in the propagation phase. Using the 



boundary conditions and the crack propagation data, the dynamic FE analysis 
determines the stress/strain field in the specimen in successive time steps. Then the 
variation of J-integral with time is obtained as shown in Fig. 4.7. 

In the beginning the crack remain stationary upto the initiation time and J- 
integral increases. At the initiation time, the crack tip starts growing under known 
dynamic displacement boundary conditions. The J -integral suddenly decreases at the 
initiation time and stabilizes after some time; the stabilized value is treated as 
propagation toughness (Jprop). 

The oscillatory behavior of J-integral is due to several reasons; 

• The free surfaces of specimen are very close to the crack tip. Some stress waves, 
emanating from the crack tip, return to the crack tip. The superposition of these 
waves and their effect on J-integral will provide ripples. 

• In the finite element modeling the crack tip is moved from one point to another 
point in discrete steps by linearly decreasing the holding back force. Thus, the FE 
code may not be simulating the crack propagation behavior very accurately. The 
nodal release mechanisms can further be modified to have smooth release of the 
crack to rule out its effect on the oscillatory nature of Jprop- 

4.3 Compilation of Results for Shorter Crack Length 

In this section, the experimental records, their subsequent analysis and the 
results from the finite element code are presented. Experimental measurements and 
subsequent finite element analysis for each experiment is shown at the end of the 
chapter. 

Two types of specimen are used in the present study. The first set employs 
specimen of a short precrack (39.5 to 42.2 mm) and cantilever thickness in the range 
of 3.5 to 4.2 mm. The second set uses specimen with precrack of longer length (49.6 
to 5 1 .0 mm); the range of cantilever thickness is same as that of the first set. 



4.3.1 Details of the Crack Lengths and Locations of Strain Gauges 

The details of crack length and strain gauge locations for different experiment 
are tabulated in Table 4.1. The thickness and width of the cantilever end are used in 
mesh generation in the FE code (Sec. 3.2.2). Taking the axes of the load bar as 
reference, the crack length and locations of velocity strain gauges are measured 
accurately with a travelling microscope. The distance between the velocity strain 
gauges is maintained between 3.0 and 4.7 mm. Also, the first velocity strain gauge is 
bonded very' close to the crack tip, about 2 mm away. 


Table 4.1 Details of Specimen and Locations of Strain Gauges for Short Crack. 


Expt 

No. 

Specimen 

Tliickness 

(mm) 

Tliickness 

of 

Cantilever 

(mm) 

Widtli 

of 

Specimen 

(mm) 

Crack 

lengtli 

(mm) 

Strain gauges locations from reference 

SGI 

(mm) 

SG2 

(mm) 

Distance 
between 
two gauges 
(mm) 

SGI 
from 
crack tip 
(mm) 

S-1 

7.6 

3.8 

24.8 

40.0 

42.7 

47.0 

4.3 

2.7 

S-2 

7.5 

3.8 

25.6 

40.1 

42.1 

45.8 

3.7 

2.0 

S-3 

7.6 

3.8 

25.1 

40.3 

42.4 

46.5 

4.1 

2.1 ^ 

1 

S-4 

8.0 

4.2 

24.5 

39.5 

41.7 

45.8 

4.1 

2.2 

S-5 

7.7 

3.9 

25.0 

40.6 

42.6 

47.1 

4.5 

2.0 

S-6 

7.8 

3.8 

25.1 

40.0 

42.8 

47.5 

4.7 

2.7 

S-1 

7.8 

4.0 

24.7 

40.4 

42.2 

45.9 

3.7 

1.8 

S-8 

7.5 

3.7 

25.5 

40.0 

42.0 

45.5 

3.5 

2.0 

S-9 

7.3 

3.8 

25.0 

40.0 

41.9 

45.4 

3.5 

1.9 

S-10 

6.0 

3.5 

25.0 

40.2 

42.6 

45.9 

3.3 

2.4 

S-11 

7.0 

3.5 

25.3 

40.1 

42.1 

45.6 

3.5 

2.0 

S-12 

7.0 

3.5 

25.0 

40.3 

42.2 

45.3 

3.0 

1.9 

S-13 

7.7 

3.5 

25.8 

42.2 

44.1 

47.1 

3.0 

1.9 


4.3.2 Crack Velocity and Initiation Time 


Table 4.2 presents the records of the velocity strain gauges responses. The 
responses of the velocity strain gauges lead to the determination of crack velocity. 
When the crack tip passes near a strain gauge the singular strain field gives a strain 




peak that is recorded in the oscilloscope (Fig. 4.4). The crack velocity is determined 
from the time difference of the peaks and the distance between the two strain gauges. 
The velocity of crack ranges between 660 to 1850 m/s. The initiation time is 
obtained by extrapolating the crack velocity upto the precrack length. The initiation 
time for the first set of experiment with short precrack length (~ 40 mm) varies 
between 44.9 and 54.9 p.s. 

TabIe-4.2 Crack Velocity and Initiation Time for Short Crack. 



4.3,3 Interlaminar Fracture Toughness 

Experimental data of the Expt. S-1 as discussed in Sec. 4.2 are summarized as 
follows : 

Voltage drop of incident pulse (Fig.4. lb) =10 mV 

Stress level in load bar corresponding to voltage drop = 259.84 MPa 



Duration of incident pulse = 12S ps 

Velocity of the load bar end (Fig.4.2) =12.5 m/s 

(between 40 ps and 120 ps) 

Maximum deflection of cantilever end (Fig.4.3) = 1.4 mm 

Initiation time (Fig.4.5) = 44.86 ps 

Crack tip velocity (Fig.4. 5) = 860 m/s 


Experimental data along with material propeities'(FL, Et, Glt, Vlt, p) are used 
to evaluate J-integral (Appendix 3). Variation of the J-integral for the stationary crack 
Expt. S-1 is shown in Figs. 4.6. The crack propagation is invoked after the initiation 
time (44.86 ps) and then the J-Integral is shown in Figs. 4.7(a) and 4.7(b). J-Integral 
for initiation and propagation is directly obtained from the figures; which are 82 and 
18 J/m^ respectively. Oscilloscope records, variation of velocity of load bar end, 
variation of deflection of cantilever end and variation of J-integral for stationary - 
propagating phase for other experiments are presented in Appendix 2. Zoomed views 
of propagating phase for all experiments are shown in Figs. 4.9 and 4.10. Crack tip 
velocity, initiation and propagation toughness are presented in Table 4.3. J-initiation 
varies from 52 to 536 J/m^ and propagation toughness 6.0 to 1 8 J/m^. 

Table 4.3 J - initiation and J-propagation Toughness for Short Crack. 


Expt. No. 

Crack Velocitv' 

a 

(m/s) 

Initiation Toughness 
J™ (J/m") 

Propagation 

Toughness 

Jprop (J/m**) 

S-l 

860 

52 

18 

S-2 

1850 

118 

7.5 

S-3 

1367 

116. 

9 

S-4 

1640 

321 

7.5 

S-5 

1286 

275 

6.0 

S-6 

1175 

418 

7.5 

S-7 

1233 

418 

11 

S-8 

1167 

384 

10 

S-9 

875 

466 

10 

S-10 

660 

414 

15 


875 

536 

15 

S-12 

750 

535 

20 

S-1.1 

667 

210 

i ^ 



4.4 Compilation of Results for Longer Crack Length 


The details of crack length and location of strain gauges for longer 
crack are tabulated in Table 4.4. The records of velocity strain gauges are tabulated in 
Table 4.5 and Jjni and Jprop are presented in Table 4.6. Jjni and Jprop were found 62 to 
400 J/m and 3.0 to 20 ]lvc^ for crack velocity ranging 633 to 2150 m/s. Zoomed 
views of propagating phases are shown in Fig.4.11 and 4.12. The oscilloscope 
records, variation of velocity of load bar end and deflection of cantilever end; 
variation of J-integral for stationary and propagating phases are presented in 
Appendix 2. 


Table 4.4 Details of Specimen and Locations of Strain Gauges for Long Crack- 


Expt 

No. 

Specimen 

Tliickness 

(mm) 

Thickness 

of 

Cantilever 

(mm) 

Width 

of 

Specimen 

(mm) 

Crack 

length 

(mm) 

Strain gauges locations from reference 

SGI 

(mm) 

SG2 

(mm) 

Distance 
between 
two gauges 
(mm) 

SGI 
from 
crack tip 
(mm) 

L-1 

8.5 

mgm 

igm 

50.2 

52.4 

56.7 

4.3 

2.2 

L-2 

7.8 

4.0 

IIIQI^III 

49.8 

52.0 

56.2 

4.2 

2.2 

L-3 

7.8 


■m 

49.6 

52.0 

56.1 

4.1 

2.1 

L-4 

8.0 

HH 

25.0 

50.6 

53.0 

56.9 

3.9 


L-5 


■■ 

25.3 

50.1 

52.4 

56.9 

4.5 


L-6 

8.3 

4.3 

25.4 

50.6 

52.8 

56.6 

3.6 


L-7 

8.4 


25.4 

50.2 

52.4 

56.6 

4.2 


L-8 

7.9 

4.1 

25.3 

51.0 

53.3 

57.1 

3.8 


L-9 

oo 

nm 

20.9 

49.7 

52.2 

56.2 

4.0 



7.4 


24.9 

50.2 

52.1 

56.6 

4.5 

1.9 

L-ll 


jjj^^ 

24.9 

50.1 

51.8 

56.0 

4.2 

1.7 

L-12 

7.1 

3.8 

25.1 

50.8 

52.4 

56.5 

4.1 

1.6 











TabIe-4.5 Crack Velocity and Initiation Time for Long Crack. 


Expt. 

No. 

Peak Responses of Strain Gauges 

Crack 
Velocity 
a = A2/AT 
(m/s) 

Initiation Time 

Ti„i = T, -A,/S 

(ps) 

SGI 

Ti (ps) 

SG2 

T: (lis) 

Time diff. 

AT “T] — T2, (ps) 

L-1 

61.5 

63.5 

2.0 

2150 

60.48 

L-2 

58.5 

62.5 

4.0 

1050 

56.40 

L-3 

62.0 

64.0 

2.0 

2050 

60.98 

L-4 

61.5 

65.5 

4.0 

975 

59.04 

L-5 

58.5 

62,5 

4.0 

1125 

56.46 

L-6 

60.0 

63.0 

3.0 

1267 

58.26 

L-7 

1 

61.0 

64.0 

3.0 

1400 

59.43 

L-8 

69.5 

75.5 

6.0 

633 

65.87 

L-9 

64.0 

66.0 

2.0 

2000 

62.75 

L-10 

61.5 

65.5 

4.0 

1125 

59.81 

L-11 

59.0 

62.0 

3.0 

1400 

57.79 

L-12 

68.0 

70.5 

2.5 

1640 

67.02 


Table 4.6 J - initiation and J-propagation Toughness for Long Crack. 


Expt. No. 

Crack Velociu 

k (m/s) 

Initiation Touglmess 

(J/m^^) 

Propagation 

Toughness 

Jprop 

L-1 

2150 

253 

10 

L-2 

1050 

125 

12.5 

L-3 

2050 

400 

10 

L-4 

975 

69 

10 

L-5 

1125 

72 

15 



L-6 

1267 

301 

10 

L-7 

1400 

242 

12 

L-8 

633 

386 

20 

L-9 

2000 

393 

10 

L-10 

1125 

38 

7 

L-11 

1400 

72 

■3, ■ 

L-12 

1640 

343 

15 







4.5 Discussion 


4.5.1 Initiation toughness of short crack 

The initiation toughness of 13 experiments with short crack is shown in Fig. 
4.13 with a best-fit linear curve. It shows that the initiation toughness decreases with 
crack tip velocity. The best-fit curve shows that Jim at crack velocity of 660 m/s is 
410 J/m^ which decreases to 200 J/m^ at crack velocity of 1850 m/s. 

Results of other investigators on glass fabric reinforced epoxy laminates donot 
exist to the best knowledge of this investigator and therefore direct comparison of 
initiation toughness under impact conditions cannot be made. However, 
Ramakhrishna (1997), Babu & Mallikharjuna (1998) and Pandey (2000) have 
conducted similar experiments on glass fibre unidirectional laminates. 

Ramakhrishna used specimen whose dimensions and precrack length were 
almost same as those of this study. Initiation toughness reported by him was varying 
between 91 and 230 J/m^. Comparison of the initiation toughness with that of the 
present study is shown in Fig. 4.14. 

While comparing these two sets of results, following differences are worth 

noting: 

• Ramakhrishna work was on unidirectional glass fibre laminates whereas this 
study is made on glass fabric reinforced polymer composites. 

• Ramakhrishna's laminates were prepared from prepregs, which were cured at 
usual temperature of 120®C. The quasistatic toughness of these laminates 
reported by him is 415 J/m^. This toughness is of the same order of laminates 
prepared by prepregs of well-known international suppliers. The prepregs used in 
the present study were made through a modified technique in which prepregs are 
cured with a cold technology using curing temperature substantially smaller than 
120®C. Consequently, the plies are bonded better with each other when a 
laminate is made. The quasistatic toughness of this laminate is experimentally 
determined to be 924 J/m^ (Appendix 3). 



• Ramakhrishna sharpened the tip of interlaminar crack of his specimen by 
extending the precrack with a wedged chisel before conducting the experiments. 
In ideal conditions, the crack front should remain normal to the axis of the 
specimen. However, some inclination from the ideal case is always present but 
efforts are made to keep the inclination minimum. It was realized, after 
conducting many experiments that in a specimen with a sharpened crack tip with 
the chisel, the probability of low inclination of the crack front is poor and 
therefore the success rate was low (<25%). Specimens used for the present study 
are as cast introduced by a thin sheet of Teflon (0.12 mm thickness). Therefore, 
the probability that crack front is normal to specimen axis is high. This leads to a 
high success rate of repeatability of experiments better than 85%. 

As expected due to low quasistatic toughness and reasons just mentioned above, 
initiation toughness reported by Ramakhrishna is small in comparison of the results 
obtained in the present study. 

The work of Babu & Mallikharjuna (1998) was similar to the work of 
Ramakhrishna. However, they used laminates of unidirectional glass fibre reinforced 
epoxy which were prepared from prepregs produced through the cold technology. 
The quasistatic interlaminar toughness was 1055 jW, which was much higher than 
quasistatic toughness (415 jW) of laminates used by Ramakrishna. Their specimen 
dimensions and precrack length were same as those of the present work. The 
specimens used by them were sharpened with a sharp nose chisel. Therefore, their 
success rate was also poor (<25%). 

Further, they used velocity strain gauges of smaller length (0.2 mm gauge 
length). Although this length of the gauges was found appropriate for experiments 
with steel (Verma, 1995) but these strain gauges are not suitable for heterogeneous 
material like polymer composites. In the present work larger strain gauge with length 
of 2 mm and width 0.9 mm and specially designed for composites are used. They are 

purchased from Tokyo Sokki Kenkyijo Co. Ltd., Japan. 

JKWTRAL LlBRAirr 

1.1. T., KAtiWJi 



Figure 4.15 shows the comparison of initiation toughness reported by Babu 
and Mallikharjuna and that of the present work. Following observations are worth 
noting: 

• Babu & Mallikharjuna conducted experiments at relatively low crack tip 
velocities in comparison of the present study. Their crack tip velocity range is 
from 472 to 933 m/s; experiments of the present work are conducted in the range 
of 667 to 1850 m/s. Average initiation toughness of the Babu and Mallikharjuna 
is 343 J/m^ and the present study is 369 J/m^ for the same range of crack velocity 
of 472 to 933 m/s. Since the results match well, crack sharpening does not seem 
to make a difference. Also the dynamic toughness of the unidirectional laminates 
and the fabric based laminates is not significantly different. The as cast crack tip 
certainly improves the reliability of success. 

• Quasistatic toughness reported by them is 1055 J/m^ is also almost same as 

* "2 

determined in the present work, which is 924 J/m . 

Pandey (2000) also worked on unidirectional glass fibre composites laminates 
made from prepregs manufacturer by the cold technology. He also used specimen as 
cast precrack introduced by a thin sheet of Teflon (0.12 mm thickness). Velocity 
strain gauges of 2 mm length and 0.9 mm width and specially designed for 
composites are used to get strain peaks in the specimen. Therefore, Pandey used 
almost same procedure for specimen preparation and velocity strain gauges that 
adopted in the present work. His results are worth comparable to the present study 
although his specimen is made of unidirectional glass fiber reinforced polymers. 

The comparison of results reported by him and the present work for initiation 
toughness of short crack is shown in Fig, 4.16. In both cases, J-initiation decreases 
with interlaminar crack velocity and the results are similar although the laminate 

material was different. 



4.5.2 J- Propagation of short crack 


The propagation toughness of 13 experiments with short crack is shown in 
Fig. 4.17 with a best-fit linear curve. It shows that the propagation toughness 
decreases with crack tip velocity. From the fit curve Jprop at crack velocity of 660 m/s 
is 14.5 J/m^ which decreases to 6 J/m^ at crack velocit}- of 1850 m/s. 

As stated in the previous subsection, the results of other investigators on glass 
fabric reinforced epoxy laminates are not available therefore direct comparison of 
propagation toughness can not be made. However, results reported by Babu & 
Mallikharjuna (1998) and Pandey (2000) for unidirectional glass fibre polymer 
laminates are used to compare the results of the present study. 

Babu and Mallikharjuna (1998) investigated with specimen geometry and 
precrack length same as of the present work. The average value of the results 
reported by Babu and Mallikharjuna is 65 J/m^ for crack velocity between from 422 
to 933 m/s. Comparison of propagation toughness reported by Babu and 
Mallikharjuna with that of the present work is shown in Fig. 4.18. The Jprop of the 
study is substantially lower than Jprop Babu & Mallikharjuna. The nature of interface 
between two like plies of unidirectional laminate is quite different from the interface 
between two plies of fabric reinforced laminate. In the unidirectional case, fibers of 
adjacent like plies lend to dig into each other. On the other hand, fibers of a fabric ply 
are well constrained. The lower Jprop of fabric laminate is probably due to this 
difference. 

Further, Pandey (2000) did work on unidirectional glass fiber composites for 
same geometry and precrack length. His specimens are also prepared through the 
cold technology. Comparison of propagation toughness reported by him and the 
present work shewn in Fig. 4.19. It is worth noting that Pandey reported low values 
of propagation toughness with average of 20.5 J/m^ whereas Babu and mallikharjuna 
reported 65 J/ml The Jprop of Panday’s investigation on unidirectional laminate is 

slightly higher than Jprop of this study on fabric laminates. 



4.5.3 Initiation and propagation toughness of long crack 

The initiation toughness of 12 experiments with long crack is shown in Fig. 
4.20. Unlike in the case of short crack where Jini decreases with the crack velocity, Jjni 
of long crack does not show any trend. The average value of the initiation toughness 
of long crack specimen is found 224.5 J/m^. The average value of the initiation 
toughness of short crack is 328 J/m^, which is higher than that of long crack. Because 
of the high scatter, the difference is not significant. 


The propagation toughness of these 12 experiments is shown fig. 4.21 with 
that of short crack. The values of Jprop and their variation with crack velocity are 
similar to the Jprop of short crack; that is, Jprop decreases with crack velocity. 

4.6 Closure 

In this Chapter, the initiation toughness of 13 experiments with short crack 
and 12 experiments of long crack was tabulated and compared with other 
investigators. These 25 experiments are conducted on fabric glass reinforced epoxy 
laminates. Specimens used for the present study are as cast introduced by a thin sheet 
of Teflon. Also velocity strain guages specially made for composites are used to get 
strain peaks. 

The variation of Jini of short crack decreases with crack velocity whereas it 
does not show any trend for long crack. However, the variation of Jprop decreases with 
crack velocity for both short and long crack lengths. The effect of crack length 
variation on Jini and Jprop is not significant. 




Fig:4.1(a) Oscilloscope records of Expt S-1 



Fig. 4. 1 (b) Zoomed view of oscilloscope records of Expt. S-1 




Fig.4,2 Velocity of load bar end obtained through incident 
-reflected pulses of the load bar for Expt. S-1 
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Fig4.3 Displacement of cantilever end obtained by integrating 
the velocity input curve for Dept. S-1 
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Fig. 4. 4 Magnified view of responses of the velocity strain gauges 
for Expt. S-1 



Fig.4,5 Extrapolation of Initiation time 
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Fig.4.7(b) Zoomed view J-pr 
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Fig.4.8 Dynamic fracture toughness variation upto initiation time 
[John and Rosakis( 1997(a)] 
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Fig.4.9 Zoomed views of J-propagation [Expts. S-2to S-7] 
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Fig.4.10 Zoomed views of J-propagation [ Expts. S-8 to S-13] 
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Fig. 4. 1 9 Comparison of J-prop with that of Pandey 
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Fig.4.20 Comparision of J-ini of long and short cracks of the study 
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Fig.4.2I Comparision of J-prop for long and short cracks of the study 




Chapter 5 

Conclusion and Scope for Further Work 


5.1 Conclusion 

The present work is carried on fiberglass fabric/epoxy laminates for finding 
dynamic interlaminar toughness under impact loading through a hybrid technique 
using experimental measurements and finite element analysis to determine the 
initiation and propagation toughness of DCB specimen. 


Experimental technique consists of a specimen; the front face of which is 
bonded to a rigid block and the other cantilever was loaded through a load bar. The 
load bar is impacted with a striker accelerated in an air gun. An incident compressive 
elastic pulse, generated in the load bar, transfers a portion of energy to the specimen; 
of the energy, the rest is reflected as the reflected pulse. The incident and reflected 
pulses are monitored by strain gauges bonded on the load bar fi'om which deflection 
of cantilever end is calculated 


The response within specimen after the impact is simulated using the FE code 
and boundary conditions determined through the experiment. After each successive 
time step, stress and displacement fields are determined in the entire specimen. An 
appropriate path is chosen to find J-integral around the crack tip. Initially the crack 
remains stationary and ] builds up. At the initiation time, the crack starts propagating; 
the crack growth is introduced in the numerical simulation by gradual release of nodal 
forces at the crack front. Further, the crack velocity was measured by monitoring the 
strain response of the two strain gauges bonded in series ahead of the crack tip. 

The initiation toughness (Jini) was found to vary between 52 to 550 J/m' for 
short crack and between 38 and 400 J/m^ for long crack. The value of the propagation 
toughness (W lies between 6 to 18 jW for shon crack and between 3 and 20 J/m^ 
for long crack. The propagation toughness was measured for crack velocity vary.ng 



in the range 660 to 2200 m/s. These results were then compared with quasistatic 
interlaminar toughness (Gic). It has been observed that dynamic interlaminar 
toughness of the crack is much smaller than the quasistatic interlaminar toughness. 


5.2 Scope of Further Work 

• More experiments are needed to study the effect of crack length over the 
wider range, 25 mm to 65 mm, on dynamic toughness. 

• The thickness of the cantilever can be varied to find its effect on the values 

of Jini and Jprop- 

• Advance digital storage oscilloscope/high speed photography may be used 
to record strain gauges responses more accurately. 
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Appendix -1 

Experimental Characterization of Composites 

Experimental characterization refers to the determination of the material 
pioperties through tests conducted on suitably designed specimens. Elastic constants 
and strengths are basic mechanical properties of materials. For a unidirectional 
laminate, there are four independent elastic constants; the elastic module in the 
longitudinal and transverse directions, the shear modulus, the major Poisson ratio. 

For a fabric composite material having equal number of yams in warp and weft 
directions, the elastic module in the longitudinal (warp) and transverse (weft) 
directions are approximately equal. Similarly, major and minor Poisson ratios are 
also equal. In the present study, the glassfibre fabric having equal number of yarns in 
warp and weft directions is used. Therefore, only three tests are needed to determine 
El ,vlt and Glt . 

DETERMINATION OF El ANDvlt 

The static unidirectional tension test is conducted to determine elastic modulus 
and Poisson ratio of the material. Straight specimen with 25.4 width and 153 mm 
long are employed with end tabs of 38 mm long on each end (Fig. A. 1.1). The 
specimen are cut from a laminate of 16 plies. 

The test are performed on Instron 850 at displacement control loading rate of 2 
mm/min. Fig. A. 1 . 1 shows dimensioning of the tensile test specimen . Strain gauges 
BELA-5-8 supplied by Toky Sokki Kankyijo Co. Ltd., Japan are used. The Poisson 
ratio is measured by placing the two strain gauges at the center of specimen in 
longitudinal and lateral directions. The responses of the strain gauges is measured 
through strain indicators. The Poisson ratio is given by 
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Table A.L1 Experimental values of El and vlt 


Expt. N 

Longitudinal Modulus 

El (GPa) 

Major Poisson Ratio 

Vlt 

1 

26.72 

0.20 

2 

26.56 

0.19 

3 

25.53 

0.21 

4 

27.79 

0.20 

5 

24.55 

0.21 


The average value of El = 26.23 GPa and average ,vlt = 0.202 were used in the 
computer simulation. 


DETERMINATION OF ,Glt 


In plane shear properties of a fabric glass composite are determined by conducting 
a tension test on a [±45]s laminate. Glt is determined for [±45]s laminate of 16 layers at 6 
mm/min displacement by using expression 


Table A.L2 Modulus of Rigidity 


^LT ~ 


2 ( 1 + 


Expt. No 

Modulus of Rigidity Glt (GPa) 

i 

3.573 

2 

3.582 

3 

3.500 

4 

3.500 






Appendix 2 

Details of Experiments 

The figures of oscilloscope records, variation of velocity of load bar end, 
displacement history of cantilever-end, variation of J-integral for stationary crack, and 
J-integral of propagating crack are shown in this appendix for Expts. No. S-2 to S-13 
and L-1 to L-12. Experimental number and concerned figures are tabulated in 
following Table; 


Expt. No. 

Concerned Figures 

S-2 

A.2. 1 (a, b, c, d, e) 

S-3 

A.2.2 (a, b, c, d, e) 

S-4 

A.2.3 (n, b, c, d, e) 

S-5 

A.2.4 (a, b, c, d, e) 

8-6 

A. 2. 5 (a, b, c, d, e) 

S-7 

A.2.6 (a, b, c, d, e) 

S-8 

A.2.7 (a, b, c, d, e) 

S-9 

A.2.8 (a, b, c, d, e) 

S-IO 

A.2.9 (a, b, c, d, e) 

S-11 

A.2.10(a,b, c,d,e) 

S-12 

A.2.11(a, b,c,d,e) 

S-13 

A.2.12{a,b,c,d,e) 

L-1 

A. 2. 13 (a, b, c, d, e) 

L-2 

A.2.14(a,b, c,d,e) 

L-3 

A.2. 15 (a, b, c, d, e) 

L-4 

A.2. 16 (a, b, c, d, e) 

L-5 

A.2.17 (a, b, c, d, e) 

L-6 

A.2.18 (a,b, c, d, e) 

L-7 

A.2.19 (a, b, c, d, e) 

L-8 

A.2.20 (a, b, c, d, e) 

L-9 

A.2.21 (a,b, c,d,e) 

L-10 

A.2,22 (a, b, c, d, e) 

L-U 

A.2.23 (a, b, c, d, e) 

L-12 

A.2.24 (a, b, c, d, e) 
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Fig. A. 2. 1(a) Oscilloscope Records of Expt S-2 
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Fig. A.2.2(a) Oscilloscope records of Expt. S-3 
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Fig. A. 2. 4(a) Oscilloscope Records Expt S-5 
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Fig.A.2.5(a) Oscilloscope records Expt S-6 
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Pig. A.2. 6(a) Oscilloscope records Expt. S-7 
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Fig. A 2.6(d) Variation of J-integral for stationary crack Expt. S-7 
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Fig.A.2.7(a) Oscilloscope Records Expt. S-8 
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Fig.A.2.7(a) Oscilloscope Records Expt. S-8 
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Fig. A.2, 7(b) Variation of velocity of load bar end Expt. S-8 
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Fig. A 2. 8(a) Oscilloscope Records Expt. S-9 
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Fig. A 2. 8(d) Variation of J-integral for stationary crack Expt.S-9 
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FigA2.9(a) Oscilloscope Records Bcpt. S-10 
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Fig. A.2. 1 0(a) Oscilloscope Records Expt. S-1 1 
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Fig. A.2. 1 1 (a) Oscilloscope Records Expt. S- 12 
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Fig. A 2. 1 1(c) Displacement of cantilever end Expt. S- 1 2 






Fig. A. 2. 1 2(a) Oscilloscope Records Expt. S- 1 3 
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Fig. A. 2. 1 3(a) Oscilloscope Records Expt. L- 1 













Fig.A.2. 14(a) Oscilloscope Records Expt L-2 
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Fig. A 2. 14(b) Velocity of cantilever end Expt. L-2 
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Fig. A2. 1 5(a) Oscilloscope Records Expt. L-3 
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Fig. A 2. 1 5(d) Variation of J-integral for stationary crack Expt. L-3 



Fig, A 2. 1 5(e) Variation of J-integral for stationary & Propagating 
crack Expt. L-3 





Fig. A. 2. 16(a) Oscilloscope Records Expt L-4 










Fig. A.2. 1 7(a) Oscilloscope Records Expt. L-5 
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Fig. A, 2, 1 7(c) Displacement of cantilever end Expt. L-5 
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Fig. A. 2. 1 19(a) Oscilloscope Records Expt. L-7 
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Fig A 2, 19(b) Velocity of load bar end Expt. L-7 
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Fig.A.2.20(a) Oscilloscope Records Expt. L-8 
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Fig. A. 2. 2! (a) Oscilloscope Records Expt. L-9 
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Fig. A.2. 22(a) Oscilloscope Records Expt. L-10 
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Fig. A.2, 23(a) OsciUoscope Records Expt L-11 


800 



V^xSty.rrVs 










Time,n s 

Fig A2.24 Oscilloscope Records Expt. L-12 
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Appendix 3 


Quasistatic Interlaminar Toughness 

1 o conijiare the dynamic interlaminar toughness with quasistatic interlaminar 
toughness experiments are conducted to determine the quasistatic interlaminar 
tiMmhness (Ciu ) through a standard test (P. Kumar (1999)). The test involves pulling a 
DCIl specimen in displacement control mode. Specimen used to determine 
c|uasistatic inleriaminar toughness is identical to that of a impact testing specimen. It 
is made of 04 laminae The load is applied until the crack extends by a small distance 
( S - 10 min) The machine is stopped until the crack becomes stationary and then the 
specinten is unkraded The specimen, now vrith a longer crack length, is loaded again. 
The compliance trf the specimen is determined through the loading curve (Fig.A.3.1). 
T he specimen is subjected to several loading - unloading compliance and critical load 
for several crack lenghs With proper data reduction scheme (P. Kumar (1999)), the 
critical energy release rate is evaluated. Value of Gic is determined from equation 


G„ 


3 A, A; 
2 B 


( 1 ) 


where Ai and Aj are given by 


ln(A,) = i[-3Iln(a) + SIn(C)] 

n 


( 2 ) 


ln(A,) = i[ln(a) + an(I’c)] 


(3) 


.-.nd P. atxl <; «c criticl load and orack length respectively. Crack length, a, 
compliance. C and critical load P, are tabulated in Table A,3.1 forExpt, 1. 



Table A.3.1 Logarithmic value of a, c and P, 


S.No 

Crack 

Displace 

Critical 

Compliance 

In(C) 

In (a) 

In (Pc) 


length 

ment at 

load 

C = U/Pc 





A 

u 

Pc 

(mm/N) 





(mm) 

(nmi) 

(N) 





\ 

B u 

2.1 ^ 

270.0 

0.00778 

-4.856 

3.433 

5.598 

7 . 

40 S 

3,1 

227.0 

0.0136 

-4.297 

3.701 

5.425 

3 

50 5 

5.(1 

212.5 

0.0235 

-3.750 

3.922 

5.358 

4 

1 60 0 

1 

6.5 

175.0 

0.0375 

-3.283 

4.094 

5.165 




Sum of column (S) 

-16.186 

15.15 

21.546 


By suing F.qs (1). (2) and (3) 

Ai 2.032 X 10'^ and 
A2 964.31 
Gk- - 944.7 J/m^ 


Similarly, one more experiment is conducted that gives 
Gk = 904,1 J/m^ 

fable A.3.2 shows the average value of Gic of the two experiments. 


Fxpt 

Specimen 

Quasistatic Interlaminar 

Average Value 

Nri 


Toughness (J/m^) 

(J/m1 

1 ^ 

{0"|m 

904.1 

924.4 

2 

101.4 

944.7 





Fig. A. 3.1 Load displacement curves. 





